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Abstract. 

We compute spectra from accretion disks around rapidly rotating 
neutron stars. The full effect of general relativity is considered for the 
structure calculation of the stars. We take into account the Doppler shift, 
gravitational redshift and light-bending effects in order to compute the 
observed spectra. To facilitate direct comparison with observations, a 
simple empirical function is presented which describes the numerically 
computed spectra well. This function can in principle be used to dis¬ 
tinguish between the Newtonian spectra and the relativistic spectra. We 
also discuss the possibility of constraining neutron star’s equation of state 
using our spectral models. 


1. Introduction 

A large number of low mass X-ray binaries (LMXB) are believed to harbor neu¬ 
tron stars, rotating rapidly due to accretion-induced angular momentum trans¬ 
fer. These systems show many complex spectral and temporal behaviors. One 
of the main purpose for studying such behaviors is to understand the properties 
of very high density (~ 10^® g cm“^) matter at the neutron star core. Such 
high densities can not be created in the laboratory and only the study of these 
sources can give a possible answer to this fundamental question of physics. Here 
we calculate the equation of state (EOS) dependent model spectra of the accre¬ 
tion disks around rapidly rotating neutron stars. These models, when htted to 
the observed spectra, can in principle help to constrain EOS models and hence 
to understand the properties of core-matter of neutron stars. 
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2. Formalism and Results 

We compute the structure of a rapidly rotating neutron star for realistic EOS 
models, gravitational masses and rotational speeds using the same procedure 
as Cook, Shapiro, & Teukolsky (1994). We use the axisymmetric metric (see 
Bhattacharyya et al. 2000 for description) 

dS'^ = + f^dO'^) + sin^ 9{d4> — codt)^. (1) 

To calculate the metric coefficients and the bulk structure parameters of the 
neutron star, we solve Einstein’s field equations and the equation of hydrostatic 
equilibrium simultaneously. Eor a thin blackbody disk, we compute the tem¬ 
perature profile (see Bhattacharyya et al. 2000) and hence the spectrum (see 
Bhattacharyya, Bhattacharya, & Thampan 2001b) considering the effect of the 
Doppler shift, gravitational redshift and light-bending effects. We calculate the 
spectra for different EOS models, which therefore gives a way to constrain EOSs 
when fitted to the observed spectra. However, the computation of the complete 
spectrum in this manner is rather time-consuming and therefore not quite suit¬ 
able for routine use. Therefore, in order to make our results available for routine 
spectral fitting work, we intend to present a series of approximate parametric fits 
to our computed spectra in a forthcoming publication. As a first step towards 
this aim, here we present an analytical function which describes both the rela¬ 
tivistic spectra and the Newtonian spectra well. This function with three free 
parameters {So,ld and Ea) is given by (see Bhattacharyya, Misra, & Thampan 
2001 a for details) 

S,{E) = SoE-^/\^yexp{-^), (2) 

where, 7 = —(2/3)(l -|- E/3/Ea), E is the energy of the photons in keV and 
Sf{E) is in units of photons/sec/cm^/keV. The value of /3-parameter is ~ 0.4 
for the Newtonian case, while it ranges from 0.1 to 0.35 for the relativistic cases 
(if the value of inclination angle is not too high). Therefore, constraining this 
parameter by fits to future observational data may indicate the effect of strong 
gravity in the observed spectrum. 

3. Conclusion 

This work is a step forward towards constraining EOS models of neutron stars 
by fitting the observed spectra of LMXBs. 
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